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ABSTRACT 

We consider the one-dimensional Hubbard model at half filling. We show 
that both excitation spectrum and S-matrix are determined by the SO (4) 
symmetry of the model. The complete set of excitations is given by the 
scattering states four elementary excitations, which form the fundamental 
representation of SO (4). We evaluate the exact S-matrix, which satisfies 
the Yang-Baxter relation. The results for the repulsive and attractive 
Hubbard model are related by an interchange of spin and charge degrees 
of freedom. 
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Strongly correlated electronic systems are currently under intense study in relation with 
high-T c superconductivity. The two-dimensional Hubbard model is in the very center of these 
investigations. It is believed to share important features with its one- dimensional analog' 1 '. 
A very interesting problem is the separation of spin and charge in the Hubbard model in 
one' 2 ' 3 ' and two dimensions' 4 '. In one spatial dimension the Hubbard model can be solved 
exactly by means of the Bethe Ansatz' 5 '. The Hubbard hamiltonian is given by the following 
expression 

H = - E E ( c U+i,. + 4 + i,^>) + ^EKt - - \) ■ c 1 ) 

j=i <r=T,l j'=i 

Here cj a are canonical fermionic creation operators on the lattice, j labels the sites of a 
chain of even length L, a labels the spin degrees of freedom, U is the coupling constant, and 
n- = cj,o- c j,°- is the number operator for spin a on site j. The hamiltonian (1) is invariant 
under a SO(4)=SU(2)xSU(2)/Z 2 algebra^ 8 ! . The two SU(2) algebras are the spin SU(2) 
with generators 

S = EsW> ^ = X>J,i<*t» S' = j2\(n jtl -n jA ), (2) 
j=i j=i j=i 

and the r]-SU(2) algebra with generators 

l l ^ L 

v = E ! { y . ^ = E' '^..o.- ' ^ 2 = 2 EKt + - 1) • (3) 

j=i j=i j=i 

All 6 generators commute with the hamiltonian (1). Below we first give a complete 
classification of the excitation spectrum of the repulsive Hubbard model by means of its 
SO (4) symmetry. The results are such that the spin and charge degrees of freedom completely 
separate from one another on the level of the elementary excitations. We then determine the 
exact two-particle S-matrix, which is of importance in the study of transport properties' 9 '. 
The S-matrix is SO (4) invariant and satisfies the Yang-Baxter relation, which implies the 
2-particle reducibility of the N-body scattering matrix. We conclude with a discussion of 
excitation spectrum and S-matrix of the attractive Hubbard model. Detailed derivations of 
the results presented in this letter can be found in [10]. 

Let us now turn to a classification of the excitation spectrum over the repulsive half-filled 
ground state (U > 0). In the many previous works on this problem (most notably [2,3,11]) the 
correct dispersions for many excitations were obtained. However, the Lie-algebraic structure 
of the excitation spectrum was not fully revealed. We classify all excitations with respect to 
representations of the SO (4) symmetry algebra. The spectrum is described by the scattering 
states of four elementary excitations only. Two of the elementary excitations carry spin 
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but no charge. This means that they transform in the spin-| representation of the spin- 
SU(2), and are a singlet with respect to the 77-SU(2). The resulting SO(4)-representation 
is denoted by (^,0). We shall call these two excitations spinons (They are similar to the 
two elementary excitations in the spin-| Heisenberg XXX antiferromagnet at zero magnetic 
fieldt 12-14 '). Energy and momentum of these particles are expressed in terms of a spectral 
parameter A. It plays the same role as the rapidity for the case of relativistic particles. 
Energy e s (A) and momentum p s (X) 



... 7r f°° duo J (u)sm(u\) ... 
2 J uo coshfwt/j 



/■oo 

e s (X) = 2 / 
Jo 



duo J 1 (w)cos(u;A) 
uo cosh(u;£/) 



(4) 



Here J n (co) are Bessel functions. This dispersion relation was obtained in [2,11]. 

The other two elementary excitations carry charge ±e but no spin. They are called holon 
and antiholon and form the (0, |) -representation of SO(4), i.e. they transform in the spin-| 
representation of the ry-SU(2) (3) and are a singlet with respect to the spin-SU(2)(2). The 
momentum of holon and antiholon differs by 7r as ry does not commute with the momentum 
operator but changes the momentum by n. The holon carries charge +e, where — e is the 
charge of the electron. Energy and momentum are parametrized by the "bare" momentum 
k h of the holon, and sin(k h ) plays the role of the spectral parameter. Energy and momentum 
are given 



r, (hh\ - n _ uh _ f°° du J (u)sm(u sm(k h ))e 
Ph{ ' ~ 2 k J uj cosh(c^) 

e h (k h ) = 2U + 2 cos(k h ) + 2 / 

Jo 



-u>U 



(5) 

du J 1 (a;)cos(ct; sin(A; /l ))e-^ v ; 

uj cosh(ujU) 



The antiholon can be obtained from the wave function of the holon by applying the operator 
. Its energy and momentum are found to be 

£ah( k ) = £h( k ) , Pah( k ) = -t + Ph( k ) ■ ( 6 ) 

The spectrum of all excitations of the repulsive half-filled Hubbard model is given by 
the scattering states of the four elementary excitations listed above. The proof of this is 
similar to the proof of the corresponding statement in the spin-| Heisenberg XXX antifer- 
romagnet [ 12_14 1. The basis of the proof is a theorem of [15], which states that (i) All Bethe 
Ansatz states with < below (or at) half-filling are lowest weight states of the SO (4) 
symmetry algebra, and (ii) The set of states obtained by acting with the SO (4) raising op- 
erators on the Bethe states is complete, i.e. it forms a basis of the electronic Hilbert space. 
This theorem implies that a complete set of excitations over the ground state is obtained 
by constructing excited states by means of the Bethe Ansatz and then act on these with 
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the SO (4) raising operators. A complete set of Bethe Ansatz excitations can be constructed 
by using the zero-temperature limit of the Thermodynamic Bethe Ansatz' 16 ^. Like in the 
XXX-antiferromagnet' 12 ' 13 ' elementary excitations always occur in pairs. This is due to 
the Z 2 -factor, which projects out all representations of SU(2) xSU(2) with half-odd integer 
quantum number: S z + r/ z = N, + |f (L is even). Accordingly there are no one-particle 
excitations, and more generally no excitations with and odd number of particles (spinons or 
holons/antiholons) in the half- filled Hubbard model. 

In order to construct the two-particle scattering matrix a complete knowledge of all 
excitations in the two-particle sector is necessary. In the two-particle sector there exist three 
kinds of excitations: Scattering states of two spinons (forming the representation (1, 0)©(0, 0) 
); Scattering states of two holons (in the representation (0, 1) © (0, 0)); and finally scattering 
states of one holon and one spinon (forming the (|, |) representation of SO (4) ). Bethe states 
in the Hubbard model are parametrized by sets of spectral parameters {A™}, {A^ 1 }, 

which are subject to the following sets of coupled algebraic equations' 16 ' 



n=l a=l ^ ' n=l a=l ^ 

N e -2M' /An • j,\ oo M m / \n _ \m 

E o - w + E E e„ 

j=l ^ ' m=l /3=l ^ 



(7) 



L [sin- 1 (A* + inU) + sin' 1 (A^ - inU 

. 2 ^/f%(^) + f;ge„,,(M) , W 

j=l ^ ' m=l (3=1 ^ ' 

Here L = 2 x odd is the even length of the lattice, J™, and J'™ are integer or half-odd 
integer numbers, N e is the total number of electrons, M' = Y^=i^'w ®( x ) = 2arctan(x), 
and 

*nm(*) = { (9) 

[ 2 0(f) + 20(f) + ... + 20(^) + 0(|;) ifn = m . 

The spectral parameters A™, A^ are the real centers of complex roots ("strings" of lenght 
n) of the Lieb-Wu equations (see [16]), M n and M' n are the numbers of A-strings of length 
n and A'-strings of length n in a specific solution of the system (7). The integer (half-odd 
integer) numbers in (7) have the ranges 



1 00 
4 n |<-(iV e -2M'-^t„ m M m -l), 



m=1 00 (10) 

1 L - 1 L - 1 

•C l<2^ - K + 2M> - J2 t nm M' m - 1) , — < I 3 ,< — , 

m=l 
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where t mn = 2min{m, n} — 5 mn . Each set of "integers" {Ij}, {Ja}i {J'a} is m one-to-one 
correspondence with a set of spectral parameters, which in turn unambiguously specifies 
one eigenstate of the hamiltonian (1). Thus the ground state and all excited states can be 
constructed by specifying sets of "integers" in (7) and then taking the thermodynamic limit, 
in which (7) reduce to coupled integral equations (see e.g. [5,2,3,11]). The half-filled ground 
state is obtained by taking M x = § and N e = L in (7) (M n = Vn > 2, M' n = Vn) and 
filling all vacancies for the "integers" /■ and J* symmetrically around zero' 5 '. This amounts 
to filling the Fermi seas of spectral parameters kj between — ir and n, and of A* between 
— oo and oo. The ground state is an SO(4) singlet (0,0). 

The scattering states of two spinons are similar to the the XXX-antiferro magnet ' 12 ~ 14 '. 
The lowest weight state of the spin-triplet (S = 1, S z = —1, rj = 0) is obtained by diminishing 
the number of down spins by one as compared to the ground state M 1 = § — 1- This 
enlarges the allowed range of "integers" J* (10) by one, so that we end up with two holes 
in the Fermi sea of the A^. We denote the corresponding spectral parameters by and 
Ag. Energy and momentum of this state are given in terms of the two spinons according to 
E ST = e s (A^) + e s (A2) and P ST = p s (A^) + p s (A%). The lowest weight theorem cited above 
guarantees that not only S z = — 1 but that also S = 1. The two other states of the triplet 
can be obtained by application of S^. 

The S = 0,i] = spin-singlet scattering state of two spinons is constructed by choosing 
M 1 — j£ — 2 and M 2 — 1. This introduces again two holes in the A^-distribution. In addition 
there is one A-string of lenght 2 present. The allowed range of "integers" for the 2-string 
(10) is only. Solving the equation (7) for the spectral parameter A^ 2 ) we find A^ 2 ) = Al ^ 2 . 
Energy and momentum of the singlet are the same as for the triplet. 

The scattering state of a spinon (with spin up) and a holon can be constructed from the 
Bethe Ansatz by choosing N e — L — 1 and M 1 = y — 1, which introduces one hole in the 
distribution of the k^s and the A^'s respectively. We denote the corresponding spectral 
parameters by k h and A h . The spin of this excitation is S — \ — 7], S z = — \ = r/ z . 
Accordingly it carries charge +e. Energy and momentum are gven by E sc = e s (A h ) + e h (k h ) 

and P sc = p s (A h ) +Ph( kh )- The other three states of the so ( 4 ) multiplet can be obtained 
by applying S\ r)\ and both. 

The scattering state of two holons with r\ = 1, rf = — 1, S = (charge triplet) is constructed 
by making two holes in the distribution of k^s (N e — L — 2). We denote the corresponding 
spectral parameters by k\ and k\. Energy and momentum of this state are given by E CT = 
e /i(^i) + e h(^2) an d Per = Ph(ki)+Ph(k2). The other two states of the multiplet are obtained 
by applying rf. 

The charge singlet S = 0,r] = = r] z scattering state is characterized by taking N e = 
L, M 1 — \ — 1 and M' 1 = 1 (insertion of a bound state of two fc's and one A'). This 
introduces again two holes in the kj Fermi sea, and in addition there is a so-called A'-string 
present. The spectral parameter of the A'-string is determined by the system (8) to be 
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A' 1 = sm ( fc i)+ sm ( fc 2) _ Energy and momentum of the ^-singlet are E cs = e ah (k^) + e^k^) and 
Pes = Pah(ki) + Ph(k%)- This completes our discription of all two-particle states. 

We now turn to the evaluation of the two-body scattering matrix. In Quantum Mechan- 
ics the scattering matrix can be extracted from the asymptotics of the wave-function of the 
scattering state' 17 l The scattering of the collective excitations described above is a quantum 
mechanical problem. However, the ground state is highly nontrivial, as are the many-body 
wave functions given by the Bethe Ansatz. In [18] it was explained how to modify the rules 
of quantum mechanical scattering theory in order to obtain the correct scattering matrix for 
Bethe Ansatz solvable models ([18] deals with the derivation of Zamolodchikov's S-matrix 
starting from the Bethe Ansatz). The same method has been applied to the case of the XXX 
Heisenberg model' 12-14 ] (see also [19]). More precisely, the method of [18] permits the eval- 
uation of the logarithmic derivative of the S-matrix, and thus the S-matrix up to an overall 
constant factor. We have fixed this factor by considering the infinite U limit of holon-holon 
scattering. In this limit the S-matrix must reduce to the one of free fermions, which is equal 
to I. This allows us to solve the problem of fixing the constant factor in the S-matrix of 
the XXX-antiferromagnet: the S-matrix for that model is given by the expression (11). The 
complete S-matrix for the Hubbard model is 16 x 16 dimensional and blockdiagonal. The 
block describing scattering of two spinons (with spectral parameters A^ and A£) is 4 x 4 
matrix of the form 



5(Af,Aj) 



r(±±^)r(i-f; 



■ 1— ifj, •• 



1 + 



r(i + f) + i ix + i 



-n 



(ii) 



where \i = X 2U 2 and 
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1/ 



(12) 



The physical region is defined as \x > (only then scattering takes place). As a function of 
the variable fi the S-matrix (11) coincides with the S-matrix of the spin-| Heisenberg anti- 
ferromagnet. This does not hold as a function of the momentum. The 4x4 block describing 
the scattering of two charged particles (holons/antiholons) with spectral parameters sin^) 
and sin^) is 



S'(sin(^),sin(^)) 



r(i^)r(i + f) 



-1 - 



-n 



(13) 



r(±±^)r(i-f) 11- 1 

where \i = — on > 0- The block describing the scattering of one spinon on one 



holon/antiholon is purely diagonal 
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S(A\sm(k h )) = -i- : ^-^(/,// = > . 14 

\ — i exp(7r^) 2U 

The last block describes the scattering of one holon/antiholon on one spinon and is given 
by the same expression (14), where now a, = sm( - k ^~ A > o. The complete S-matrix satisfies 
the Yang-Baxter equation, which implies that the many-body scattering matrix factorizes 
into the product of two-body S-matrices. Clearly the S-matrix reflects the SO (4) symmetry 
of the model. The analytical structure of the S-matrix is as follows: The allowed strip for 
rapidities is |Im(A /l )| < U, and |Im(A; /l )| < 2U. The S-matrix has poles in the complex 
rapidity plane, but all of them are unphysical and do not correspond to bound states. 

The structure of elementary excitations in the attractive Hubbard model is as follows: 
there is a total of four elementary excitations, grouped into doublets of the spin and r7-SU(2) 
respectively. The chargeless spin-carriers have the dispersion 



p s {k) = n -k- 

Jo 



duj J (u;)sin(a; sin(/c))e 



uj cosYi(ojU) , . 

m n . TT . „ s n f°° duj JJu)cos(u sin(k))e-^ ( " o) 

e s (k) = 2\U\+2cos(k) + 2 ^ v .V , 

Jo lu cos\i(luU) 

which up to the shift of momentum by | is the same as for the holon in the repulsive case. 
The spinless charge-carriers (one particle, one hole) have the dispersions 

f°° du> J (a;)sin(a;A) h 
pP(\) = n- ' \ = n — Pc (A) , 

(16) 



o 



J oo cosh(a;C/) 



which are (up to a constant in the momenta) the same as the spinon dispersions in the 
repulsive case. The S-matrix in the attractive case coincides with the one for the repulsive 
case £ls Si function of the unformising parameter ji. There also exists a naive argument for 
the relation of attractive and repulsive Hubbard model by an interchange of charge and spin 
degrees of freedom. It is a well known fact that attractive and repulsive Hubbard model are 
connected by the particle-hole transformation (for spin-up only) 

c ,,t ^ (-i)Ht > ( 17 ) 

under which the hamiltonian (1) transforms according to H(U) — > H(—U). The transfor- 
mation (17) effectively interchanges spin and charge degrees of freedom, e.g. transforms the 
spin-SU(2) (2) into the r/-SU(2) (3) and vice versa. 

In this letter we classified the excitation spectrum and the Scattering matrix of the 
one-dimensional Hubbard model by means of its 50(4) symmetry. This emphasizes the 
importance of symmetries in systems of strongly correlated electrons. The S-matrix is an 
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important tool in the study of transport properties^. The classification of elementary ex- 
citations into spinons (which carry spin but no charge) and holons/antiholons (which carry 
charge but no spin) gives a precise meaning to spin-and charge separation in one and two- 
dimensional models of correlated electrons. In integrable models the S-matrix is often pro- 
portional to the R-matrix. In [21] a R-matrix for an embedding of the Hubbard hamiltonian 
in a commuting family of transfer matrices was derived. This R-matrix, which belongs to 
a fundamental model (where L-operator and R-matrix are essentially identical' 19 '), is differ- 
ent from our S-matrix. It would be interesting to search for a different embedding of the 
Hubbard hamiltonian, using a SO (4) invariant R-matrix proportional to our S-matrix. The 
L-operator for such a construction would have to be different from the R-matrix. 
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